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We present a quasi-particle model that allows for a simple description of the electronic properties of 
metallic helimagnets. For weak quenched disorder, we find a leading linear temperature dependence 
of the electrical conductivity for 3-d materials. This is reminiscent of the behavior of nonmagnetic 
2-d systems, and reflects a general tendency of certain properties of bulk helimagnets to appear 
effectively 2-d. The sign of the effect is opposite to that in nonmagnetic 2-d materials. These 
surprising predictions should be observable in weak helimagnets. 
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Helimagnets (e.g., MnSi, or FeGe) are magnetic mate- 
rials where the magnetization shows ferromagnetic order 
in any plane perpendicular to a certain axis, but points in 
different directions depending on the position along that 
axis, forming a spiral. The metallic hclimagnct MnSi is a 
well-studied material with unusual properties. It shows 
helical order below a critical temperature T c rs 30 K with 
a helix wavelength 2n/q w 180 A [![. Hydrostatic pres- 
sure decreases T c until the long-range order 0] disappears 
at a critical pressure p c f=a 14kbar [3J]. At higher pressure, 
there is a phase or region where short-range helical or- 
der persists and the electrical resistivity p shows a T 3 / 2 
temperature dependence in a range between a few mK 
and several K [4] . If this were the true asymptotic low-T 
behavior, it would represent non-Fermi liquid behavior 
in a bulk material, which would be very remarkable. 

The samples used in Ref. have a residual resistivity 
of about 0.3/iSlcm, which corresponds to a rather large 
elastic mean-free time r. At the experimentally attain- 
able temperatures, this places them in the weak-disorder 
regime, Tt 3> 1 [5], where the transport is governed by 
free-electron motion in between rare scattering events, 
as opposed to diffusive motion in the opposite limit, 
Tt <C 1. The weak-disorder regime has been investigated 
by Zala et al. @ for the case of electrons interacting via a 
screened Coulomb interaction in non-magnetic 2-d met- 
als, where they found a linear T dependence of the re- 
sistivity. (In 3-<i,the corresponding resistivity correction 
is of 0(T 2 InT) [7(.) Although the weak-disorder regime 
is never the true asymptotic low-temperature regime, it 
thus can display very unusual behavior. Furthermore, de- 
pending on the impurity concentration, it can represent 
the low-T asymptotics for practical purposes. 

In a helimagnet, massless fluctuations of the helix (heli- 
magnons) that couple to the electrical conductivity are an 
obvious possible source for unusual transport behavior, 
and it is natural to first study their effects in the ordered 
phase, before trying to understand their potential rami- 



fications more generally [8j. It was recently shown that 
helimagnons do indeed lead to a non-analytic tempera- 
ture dependence of the resistivity in the ordered phase 
0. However, for the most interesting observables in a 
clean system these effects provide corrections to the usual 
Fermi-liquid behavior: the leading hclimagnon contribu- 
tion to the resistivity has a T 5 / 2 behavior, and the spe- 
cific heat goes as T 2 , although the single-particle relax- 
ation rate goes as T 3 / 2 . 

In this Rapid Communication we investigate the re- 
sistivity of 3-d helimagnets in the weak-disorder regime 
(which is defined slightly differently than in the Coulomb 
case, see below). We show that the leading tempera- 
ture correction Sp(T) to the resistivity is proportional to 
T due to scattering by helimagnons. Remarkably, this 
is much stronger than either the T 5//2 clean helimagnon 
contribution or the T 2 Fcrmi-liquid contribution, and it 
is reminiscent of the behavior of 2-d nonmagnetic met- 
als. In contrast to the latter, however, the sign of the 
T-dependence is antilocalizing: there is a T-independcnt 
part that decreases the resistivity, and the T-dependent 
part is Sp oc +Tr. 

Transport theory for helimagnets is rather compli- 
cated, even in the clean case, if the Kubo formula is 
evaluated using the usual plane-wave basis. In this basis, 
the electronic Green function is not diagonal in either 
wave vector space or spin space, which makes for very 
cumbersome calculations [9| . Including impurity scatter- 
ing in this formalism would be hard. However, several 
features of this theory, which we quote below, suggest a 
much simpler effective description. The poles of the elec- 
tronic Green function, or quasi-particle (QP) energies, 
are given by 
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Here = e& — ep, with ek the electronic energy- 
momentum relation and ef the chemical potential. For a 



cubic crystal, such as MnSi, 
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with TO e the effective electron mass, and v = 0(1) a di- 
mensionless measure of deviations from nearly free elec- 
trons. We choose units such that H = ks = e = 1. 
q = (0, 0, q) is the helix pitch vector, which we choose 
to point in the ^-direction, and A is the Stoner splitting. 
The two signs of the square root represent the two Stoner 
bands, and for q = one recovers the usual ferromagnetic 
result in Stoner approximation. Excitations between the 
two Stoner bands will be gapped by A, and we can thus 
restrict ourselves to a single band by considering spinless 
QPs with a resonance frequency u>i(k). 

For the electronic QP Green function one thus expects 



(3) 



where p = (p, iuj n ) with uj n a fermionic Matsubara fre- 
quency [l(| ■ We have shown that such a QP Green func- 
tion can indeed be derived from a canonical transforma- 
tion to suitable fermionic degrees of freedom rj(p) [111 ]. 
The soft, or massless, helical degrees of freedom, the he- 
limagnons, are fluctuations of a generalized helical phase 
<j) with a frequency 



w (p) 



CzP' z 
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Here p z and p± are the longitudinal and transverse com- 
ponents, respectively, of the wave vector with respect to 
the pitch wave vector q, and the elastic constants c z and 
cj_ are given by 



7zA 2 q 2 /k 



c_l = 7±A 2 //Cf 



(4b) 



7 2 and 7j_ are model-dependent numbers. For the model 
considered in Ref. 0], their values are "/ z — 1/36 and 
7_L = 1/96. Equation (l4aj) is valid for \p\ < q. Notice 
that the helimagnon dispersion relation is anisotropic: 
it is ferromagnet-like in the longitudinal direction, but 
antiferromagnet-like in the transverse direction. The he- 
limagnon susceptibility is given by 



X (k) = (ct>(k)<t>(-k)) = 



1 



q 2 /Skl 



2N F w 2 (fc)- (in n ) 2 



(5) 



with Q n a bosonic Matsubara frequency, k = (fc,if2 n ), 
and Nf the density of states at the Fermi level. We see 
that frequency or temperature scale with the soft wave 
vector k as T ~ k z ~ k\. 

Since <f> is a phase, only the gradient of cf> is of physical 
significance and will couple to the QPs. Furthermore, the 
results of Ref. show that the leading coupling is not to 
the QP density, but rather to the z-_L components of the 
QP stress. Upon integrating out 4>, one thus expects an 
effective potential for the interaction of QPs by means of 
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FIG. 1: Effective interaction between QPs (solid lines) by 
exchange of a helimagnon (dotted line). 



exchange of hclimagnons that is given pictorially in Fig. 
[TJ and analytically by 

V oS = -A 2 ( g 2 /8m 2 ) X (fc)7(fc,p) 7 (fc,p'), (6a) 

with a vertex function 

j(k,p) = u Pz { Pl _ •fc i )/Afc 2 . (6b) 

Finally, we expect static impurities or quenched dis- 
order to couple to the QP density. Putting everything 
together, we now have the following effective action: 



S — Sq + Sint + Sd 



(7a) 



where 



So = (T/V) J2 V(p) k - (P)\ V(P) (7b) 



describes free QPs, 

k p,p' iuj 

xr](p + k,iu> + id) fj(p', iuf) Tj{p' — k,iu — ift) (7c) 



describes the effective QP interaction via the exchange 
of a helimagnon, and 

Sdis- (-l/V 2 )J2u(k-p)Tj2v(k,iLu) V (p,^) (7d) 

describes the quenched disorder. Here u is a delta- 
correlated point-like random potential with a Gaussian 
distribution whose second moment is given by l/2nN-pT. 

The Eqs. J7]) represent an effective model for QPs in- 
teracting via a dynamical potential that allows for an 
evaluation of the Kubo formula in complete analogy to 
electrons interacting via a dynamically screened Coulomb 
interaction. The only unusual aspects are, (1) the cou- 
pling to the QP stress rather than the density, and (2) 
the anisotropic nature of the exchanged excitations. The 
former is embodied in the vertex functions 7 in Eqs. ([6]) 
and easy to handle technically, but has important physi- 
cal consequences. The latter makes the bulk system be- 
have in some respects like a 2-d system. We have derived 
this effective model from the one studied in Ref. by 
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means of a canonical transformation ll| ; the above con- 
siderations represent just plausibility arguments. 

Transport theory now proceeds in analogy to the 
Coulomb problem. In order to calculate the static electri- 
cal conductivity tensor a^ we evaluate the Kubo formula 



l/T 

limn^oRe-^- / dr e m " r (T T j t (k, T ) 



in 



n Jo 



x Jj(k,T = 0) 



Iifi„^n+i0, fe=0 ' 



(8) 



with j(fe,r) the current operator in imaginary-time rep- 
resentation and T T the imaginary-time ordering operator. 
The tensor <t^ is diagonal and has two independent ele- 
ments, a L = c zz and a± = a xx = a yy . 

In the clean limit, r — > oo, we reproduce the results of 
Ref. S an infinite ladder resummation of both the self- 
energy and vertex-correction contributions to the con- 
ductivity shows that the leading low-T terms cancel be- 
tween the two sets of diagrams. The transport relaxation 
rate, and hence the resistivity, shows a T 5 / 2 behavior, 
whereas the single-particle relaxation rate goes as T 3 / 2 . 

In the weak-disorder limit, it is convenient to include 
the disorder in the Green function in the Born approxi- 
mation, i.e., to use a Green function 



G(p) = 1/ (iuj„ - oji (p) + i sgn (w„)/2r) 



(9) 



instead of Eq. ((3|) . The diagrams that contribute to the 
correction to the conductivity to leading order in the dis- 
order in the weak-disorder limit are shown in Fig. [2 they 
are the same as in Ref. 0. 

Diagrams (i) and (ii) also contribute in the clean 
limit, where their respective ladder resummations must 




FIG. 2: Leading contributions to the conductivity in the 
weak-disorder limit. Solid lines represent the Green function 
given by Eq. Q, dotted lines represent the effective poten- 
tial, Eqs. ([6]), and dashed lines represent the impurity factor 
1/2ttNft. Obvious symmetric counterparts of all diagrams 
except (ii) are not shown. See the text for additional infor- 
mation. 



be taken into account. In the presence of disorder, their 
leading terms are proportional to r 2 T 3 / 2 . (This is the 
result of adding the clean single-particle rate, which is 
proportional to T 3 / 2 , to 1/r according to Matthiessen's 
rule, and expanding.) These terms cancel between dia- 
grams (i) and (ii), as they must according to Ref. 0. The 
calculation shows that the small parameter for the disor- 
der expansion is e = l/^r 2 Te 2 /A. This is different from 
the Coulomb case, where the small parameter is 1/tT 
@. Consequently, the next-leading terms are of order 
tT . These also cancel between the two diagrams. 

For the remaining diagrams it suffices to calculate their 
leading behavior, and they can be classified with respect 
to their momentum structure. They each contain six 
Green functions that factorize into two correlation func- 
tions formed by n and 6— n Green functions, respectively, 
with n = 3 or n = 4. Power counting shows that the 
(4, 2) partitions are smaller than the (3, 3) partitions by 
a factor of (q/fcp) 2 - Dropping the former leaves us with 
diagrams (iii), (iv), (vi). Diagram (x) also is a (3,3) par- 
tition, but its leading contribution vanishes due to reality 
requirements. Furthermore, power counting reveals that 
for erj_, only diagram (iii) is of 0{tT), whereas for cl one 
needs to consider diagrams (iv) and (vi) as well. 

We thus conclude that the leading contributions to 
both the longitudinal and the transverse conductivity 
corrections are of 0{tT). We find 



<5<7l = 3<5ct i = 



-o"o 



192^/7-7: 



T 



(10) 

with do = (2fcp/37r 2 ) epT the Drude conductivity. In 
writing this result we have dropped a T-independent UV- 
cutoff dependent contribution to 5a whose sign is oppo- 
site that of the T-dependent contribution. 

We now discuss this result. First, notice the sign of 
the result. 8a increases with decreasing T, which means 
the effect is antilocalizing. This is in contrast to the 
Coulomb case, where the sign is localizing. The sign is 
the same for both the self-energy diagram (iii), and the 
vertex correction diagrams (iv) and (vi). In the former, 
it is a result of the angular dependence of the effective 
interaction, Eq. (|6b[) . which in turn is a result of the 
stress coupling, and which suppresses backscattering in 
diagram (iii). For a density coupling, diagram (iii) would 
be localizing. 

Second, we discuss the range of validity of our result. 
The weak-disorder temperature regime is bounded below 
by the requirement that the expansion parameter e be 
small, i.e., Tr 2 e|/A ^> 1. This defines a temperature 
scale T* = A/(epr) 2 . It is bounded above by the weak- 
disorder correction to the relaxation rate crossing over 
to either the clean-limit rate l/r c i can oc T 5 / 2 [91, or to 
the Fermi-liquid rate 1/tfl ~ T 2 /eF- For the former, 
we have a crossover temperature Ti_ 5 / 2 = A/(£ft) 2 / 3 , 
for the latter, Ti_ 2 = (l/T)( 9 /fc F ) 3 (e F /A) 2 . The weak- 
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disorder rate thus yields the dominant T dependence of 
the resistivity in the regime 

T*«T«Min(T 1 _ 2 ,T 1 _ 5/2 ), (11) 

and T 1 _ 2 /T 1 _ 5/2 = (q/k F f(e F /X) 3 /(eyr) 1 / 3 . If A < e F , 
then Ti_ 2 T 1 _ 5 y 2 . In a weak helimagnet, where A/e F 
might be as small as q/kp, we still have Ti_ 2 < T\_ 5 / 2 on 
account of the factor l/(e F r) 1 / 3 . For realistic parameter 
values, the weak-disorder regime will thus be 

T* « T « Ti_ 2 = T* (e F r)(q/fc F ) 3 (e F /A) 3 . (12) 

The upper limit depends strongly on e F /A, which is often 
not well known. For weak helimagnets, with A <C e F , it 
is possible to have 3> T* by virtue of e F T 3> 1. For 
larger values of A the weak-disorder correction may never 
dominate unless the system is extremely clean. In that 
case, one will have to subtract the Fermi-liquid T 2 be- 
havior in order to observe the weak-disorder contribution. 
For MnSi, with e F w 23, 000 K and e F r w 1,000 for the 
samples of Ref.|j, one has T* » 10 mK. q/k F w 0.02, but 
the value of A is not well known. For the upper and lower 
limits used in Ref. 0, namely, A = e F /2 and A = 540 K, 
respectively, we find T x _ 2 w 0.06 T* and T x _ 2 w 500 T*, 
which again underscores the strong dependence on e F /A. 
Generally speaking, weak helimagnets are the best candi- 
dates for observing the weak-disorder transport regime. 

Third, the anisotropic helimagnon dispersion relation, 
Eq. (|4ap , leads to a 3-d helimagnetic system being quali- 
tatively similar to a 2-d one with just a Coulomb interac- 
tion. The components of the soft wave vector k scale as 
k z ~ k\ ~ T. With f(k) a generic function, the relevant 
integrals that determine observables are of the form 

J dk z J dk ± k\ S(fl 2 ~k\- k\) f(k z ,k x ) oc 

f dk ± k\ e( " 2 = k ] ] f(k z = 0, k x ), 
J x ^n 2 - k 4 ± v ' 

and the dependence of / on k z can be dropped since it 
does not contribute to the leading temperature scaling. 
The prefactor of the k dependence of / is of O(l) in a 
scaling sense. As a result, the 3-d fc-integral behaves ef- 
fectively like the integral in the 2-d Coulomb case [l2j ]. 
In addition to the conductivity, these considerations ap- 
ply to weak-disorder corrections to the spin susceptibility, 
and to the specific heat [ill ]. 

Finally, we point out a generic feature relevant for the 
explanation of the T 3 / 2 dependence of the resistivity ob- 
served in the disordered phase of MnSi. The prefactor of 
T/e F in Sa^/ao, Eq. (|10p . is very small. Assuming v = 1 
and the above parameter values for MnSi, it is about 
3 x 10~ 5 for A = e F /2, and about 10~ 2 for A = 540 K. 
The small prefactor reflects the long wavelength of the 
helix on a microscopic scale. The prefactor of the T 5//2 
behavior in the clean limit is also small, for the same 



reason. By contrast, the prefactor of the experimentally 
observed (T/e F ) 3 / 2 behavior of the resistivity in the dis- 
ordered phase is of O(10 6 ). Therefore, the processes caus- 
ing the T 3 / 2 behavior in MnSi must occur on short length 
scales, and are unlikely to be related to the helical order. 

In conclusion, we predict that in the ordered phase of 
metallic helimagnets, at low temperatures in the weak- 
disorder regime, the leading temperature correction to 
the Drudc conductivity is proportional to TV, and the 
size of this effect depends on the direction of the current. 
This small but very remarkable effect should be most 
easily observable in weak helimagnets. 
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